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Abstract 

I demonstrate that the short distance contribution to K — > 7T7t decays must be 
supplemented with large distance effects. A hybrid calculation is outlined based on 
QCD diagrams supplemented by chiral contributions and 7r-7r phaseshifts. 
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At this meeting the KTeV group reported a new and smaller value of e' / e. The experimental 
values of the two groups are now consistent with each other 

e'/e = (15.3 ± 2.6) x 10" 4 NA48 [1] (1) 

e'/e = (20.7 ± 2.8) x 1(T 4 KTeV § . (2) 

On the theoretical side there were, for a long time, two types of predictions: 

i) those which predicted small, almost unmeasurable, or negative values, and 

ii) those which predicted large, positive and observable values. 

It is interesting to examine the reasons that led to the diverse predictions, and also the 
changes required in order to bring them into agreement with the experiment. For my part, 
since I belong to the second group, I will present our efforts and work that led to large 
values for the ratio. The starting point of the analysis is the effective AS = 1 Hamiltonian 
derived from QCD. The Wilson coefficients in the Hamiltonian have been calculated to 
NLO and I shall use them. It has become evident that tree diagram values for the matrix 
elements of operators are not sufficient and we need to include loop corrections, which are 
in fact unitarity corrections. The study and calculation of the loops modify the values for 
the amplitudes and improve the accuracy of the predictions. All groups which included 
loop corrections to matrix elements predicted large and positive values for {e'/e). Now 
that this fact has been recognized we face the next problem of estimating the uncertainties 
and improving the reliability of the results. This is a difficult problem because it requires 
defining a low energy limit of QCD, including loops. 

The first article to include loop corrections to matrix elements |§ studied the operators 
Qi and Q2 and found an enhancement of the 1 = amplitude. Subsequently, it was 
recognized that loops are also important for the operators Qq and Q$ [|J, because they 
enhance (Qe)o an d reduce (Qs)2 0, thus eliminating the cancellation between the two 
terms and predicting a large ratio {e'/e). Loop corrections have also been calculated in 
the chiral quark model |J with similar results. Finally, they were calculated in chiral 
perturbation theory, paying special attention to the matching of the chiral loops to the 
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QCD-scale residing in the Wilson coefficients. 



The fact that loop corrections are important is well accepted now. The remaining issue 
is the accuracy of such calculations, especially the dependence of the results on the cut-off 
of the chiral theory, since chiral theory is not renormalizable. It was shown in article |§ 
that to 0(p 2 ) the infinities of the factorizable diagrams are absorbed into the renormalized 
coupling constants 

L\ = 2.07 x KT 3 and L r 8 = 1.09 x 10~ 3 . 

This happens for all operators - those that include currents, as well as those with densities 
of quarks. The non-factorizable terms bring to 0(p 2 ) a quadratic dependence on the cut- 
off, which makes the predictions sensitive functions of the cut-off. In the meanwhile, the 
experience gained in the calculations, together with the fact that we must account for five 
physical quantities |Ao(mjc)|, |v4 2 (m^-)|, #o(m#), ^(m^-) and (e'/e)^-, allow me to abstract 
four general properties that must be included in explanations of the quantities. The aim 
here is to arrive at a low energy realization of QCD which reproduces current algebra at 
the 1-loop level, since current algebra is already satisfied at the tree level of chiral theory. 

I summarize the properties as four benchmarks. 
Benchmark 1 . The imaginary parts of the amplitude can be calculated by unitarity to be 



ImA(K+ - tt+tt ) = \A(K+ - 7rV)| (-]-) ^-^/l-^? (3) 



4ttF2 V " m 2 K 



and similarly 



„n M K = \ M K m)\ \ <-° F ^ f^- (4) 



where |^4( J FC + — > 7r + 7r )| and |Aq(X — > 7r7r)| are absolute values of the complete decay 
amplitudes for K + — > 7T + 7r° and the 1 = decay, respectively. Taking the ratios ImA/\A\, 



one arrives at 



„ n 1 m\ — 0.5m 2 / Ami , , 

sin5 o = 7 — A F2 ~ \l ^ = 0.46 5 



,9 1 ml- — 2m 2 / 4rnl , ^ 
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which give the following values for the phaseshifts 

5q = 27.4° to be compared with 34.2° from experiments, and 

8$ = —11.5° to be compared with —7° to —12° from experiments. 

One must conclude that for s = m 2 K the absolute values of K — > tttt amplitudes, combined 
with the phase space integration (unitarity), gives acceptable phases. It still remains to 
compute the real parts of the amplitudes at s = m^, as well as real parts of the matrix 
elements for various operators. 

Benchmark 2 . Practically all articles use the effective AS = 1 Hamiltonian derived from 
QCD. At a high momentum scale, like // = 2 GeV, the Hamiltonian should give reliable 
estimates which are then extrapolated to lower momenta. In addition, at the high //-scale, 
the Wilson coefficients in the various schemes are close to each other 0. Since the quarks 
carry large momenta, I will use the lowest order contribution of quark operators between 
hadrons. The numerical estimates of such calculations are shown in Table 1, where I give 
numerical values for two momenta - // = 2.0 and // = 1.0 GeV. The change between the 
two scales is relatively small, so I have also included in the table the experimental values for 
comparison. The theoretical estimate for Ao(fi) gives only 30% of the experimental value 
and for A 2 (fi) it gives 150% of the experimental value. Thus additional (non-perturbative) 
corrections must be large and they must increase A (mK) and decrease ^(m^). The 
corrections found by the Dortmund [|], [| and Trieste H] groups, in chiral theory, increase 
(Qe)o an d decrease (Q%)2 in agreement with the above requirement. Similarly, the n-n 
phase shift for / = is positive and for / = 2 negative, bringing corrections with the 
required sign |8|, ||, [TIJ . 





Ao(fj) in GeV 


A 2 {n) in GeV 


2 GeV 


0.62 x 10" 7 


0.25 x lO" 7 


1 GeV 


0.83 x 10~ 7 


0.23 x 10~ 7 


Exper. 


3.33 x 10~ 7 


0.15 x 10~ 7 



Table 1: Contribution from QCD alone 
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Benchmark 3 . It is evident from Table 1 that the change of the scale /i is not large 
enough to provide the required corrections. In fact, several operators of QCD have a 
similar development between 2.0 and 1.0 GeV. The running of C§{q 2 ) (dashed curve) and 
C§(q 2 ) / m 2 s (q 2 ) (solid curve) in this momentum region are plotted in figure 1. The variation 
of the ratio is much smaller than the running of C§(q 2 ) alone. 
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Figure 1: Evolution of Ce(q 2 ) (dashed curve) and C§(q 2 ) / m 2 (q 2 ) (solid curve) as functions of q 



This stability was already noticed in |TT| and a similar one appears for C%{q)jm s {q) 



Operators whose anomalous dimensions are zero are called marginal and they remain 
constant over extended regions of momenta. A similar property appears, numerically, in 
the sum of coefficients [|7| 

Zi(/i) + Z 2 {fj) ~ 0.7 to 0.8 for 0.6 < /x < 2.0 GeV . (7) 



These results from QCD suggest that marginal quantities of QCD could be extrapolated 
to low momenta. In this case, large corrections must reside on the values of the matrix 
elements as non-perturbative effects. 
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Figure 2a: The phaseshift 8q as function of 
the energy E = ^/s 
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Figure 2b: Inelasticity t]q for the 1 = 
channel as a function of the energy 



Benchmark 4 . We have a lot of data on n-n scattering, some of them being shown in figures 
2a and 2b. The data indicate that the phaseshift for the 1 = channel remains elastic 
up to 900 MeV and then the resonance /b(980) appears, which forces the 5$ phaseshift to 
pass through 180°, as shown in figures 2a and 2b []I2"| . The behavior of 5q must help the 
convergence of the 1 = amplitudes, since sin 5° goes through zero at E = 1 GeV. The 
1 = 2 phaseshift is negative, smaller and has no special structure. The dynamics of the 
7r-7r scattering play a role in the decay amplitudes and must be included in the analysis. 
Contributions from the phaseshifts will bring corrections with the desired signs: increase 



(Qe)o an d decrease {Qs}2- Suggestions along these lines have already been discussed [|| [TO 



It is suggestive that an explanation of the ^-decays must include the above properties. 
There are interpolations suggested by the above benchmarks. I will outline a method in this 
context, which tries to improve the high energy region of chiral calculations (p ~ lGeV). 
I will replace the i-T-meson with the divergence of the AS = 1 axial current and try to 
apply current algebra. The divergence of the axial current is realized below 1.0 GeV by its 
chiral representation and above that energy by the quark representation. The amplitude 
is an analytic function of the invariant mass of the two-pions-squared, s, and satisfies a 
dispersion relation 



_ , , , 1 f g o Disc. A^s) , 1 r m % Disc. AAs) , , o , 
ReA I (a) = - ds + - — ds . (8) 

7T Jim? S — a 7T J s S — O 

The value of sq is an intermediate scale, with y'so = 1-0 to 1.4 GeV. The splitting of the 



integral is intentional so that for the high energy term I use the QCD amplitude. For 
instance, consider the amplitude 

A Q(i (s,fi) = C 6 (s^)(irrr\Q 6 \K} (9) 

with the discontinuity given by [|13|, 14] 



Disc.A Q6 ( S ,/i) = I -^L— - CM (intjpj {im\Q,\K). (10) 

In this energy region the momenta in Qq are large and the tree contribution for the matrix 
element will suffice. The integral over the discontinuity is finite and no subtraction is 
required |14 . 



The discontinuity in the low energy integral is assumed to be |J 

Disc.A Q6 (s) = -4v/3L 5 (^^) ^|go(s) (11} 



with the Goldberger-Treiman factor |L5 



T] sin 25° cos 5° 

9o{s) = —— , 12 

I + 770 cos 2d u 

r)o the inelasticity and 5° the phaseshift. This representation is certainly valid below the 
inelastic threshold. Since the 1 = channel remains elastic up to 1.0 GeV, I will use this 
form. The enhancement factor from the first integral in eq. (|j) is given by 

/ \ 1 f*° 9o(s) m^(lGeV) ,„ , „ , , 10 . 
F 6 (s ,m k ) = - ^ I (l±h(s)) ds. (13) 

7T Jim 2 S — 171% mi(s) 



The phenomenological function h(s) = | y ~i'^ 5 J m the integrand is introduced in order 
to study the sensitivity of the enhancement factor on off-the-mass shell effects [[L6| (see 
discussion). Figure 3 shows the enhancement factor as a function of E = t/so. The factor 
F 6 (s , m>k) increases up to 900 MeV and then flattens out. To obtain the solid curve I set 
h(s) = and the strange quark mass in the denominator equal to m s (l GeV). The dashed 
and dashed-dotted curves include the running strange quark mass and ±h(s), respectively. 
The enhancement factor lies between 0.40 and 0.55. A similar calculation for Fs(so,mk) 
gives a depletion factor with the approximate value of —0.20 to —0.30. 
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Figure 3: The enhancement factor ^(so, m-fc) as a function of E = y/so 

The loop corrections increase (Qe)o by 50% and decrease (Qs)2 by 20% to 30%. For 
these values the ratio attains values 

(e'/e) « (15... 20) x 10" 4 , (14) 

which do not include the isospin breaking term yet. The range given in eq. (|nj) includes a 
rough estimate of the theoretical uncertainties. A detailed calculation with a description of 
the current algebra and a precise estimate of the theoretical uncertainties will be presented 
in a publication. 
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